Formation of Pairing Fields in Resonantly Coupled Atomic and Molecular 

Bose-Einstein Condensates 
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In this paper, we show that pair-correlations may play an important role in the quantum statis- 
tical properties of a Bose-Einstein condensed gas composed of an atomic field resonantly coupled 
with a corresponding field of molecular dimers. Specifically, pair-correlations in this system can 
dramatically modify the coherent and incoherent transfer between the atomic and molecular fields. 
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In quantum field theory, superfluidity and long-range 
order are usually identified with a complex-valued order 
parameter defined as the expectation value of the field 
operator. In the case of a dilute Bose gas, the evolu- 
tion of the order parameter is given in the mean-field 
approximation by the well-known Gross-Pitaevskii equa- 
tion. Recently, this equation has been extensively applied 
to studies of the zero-temperature collective behavior of 
metal-alkali Bose-Einstein condensates confined in mag- 
netic traps 0. Quantitative agreement has been found 
between theory and experiment on a wide range of dis- 
tinct phenomena (e.g. Refs. 

Although identical in form to a non-linear Schrodinger 
equation, it is notable that the Gross-Pitaevskii equation 
does not by itself describe the full evolution of the quan- 
tum state. In general, a complete description requires 
knowledge of the evolution of a hierarchy of correlation 
functions, and not just the evolution of the expectation 
value of the field. In particular, in low temperature sys- 
tems, it is well known that pairing fields associated with 
two-particle correlations can play an important role. In a 
quantum degenerate Fermi gas, pairing may radically al- 
ter the equilibrium properties, giving rise to superfluidity 
in 3 He and superconductivity in electron systems. Corre- 
lations can also be important in Bose systems. For exam- 
ple, squeezed states of light — formed when photons are 
generated in pairs in nonlinear media — have been studied 
extensively in the subject of quantum optics 

In this paper, we investigate the role of pair- 
correlations on the macroscopic dynamics of a di- 
lute Bose-Einstein condensate composed of atoms and 
molecules which are resonantly coupled p-Hh Such a cou- 
pling may be generated experimentally p|JlO| by tuning 
the strength of an external magnetic field in the proxim- 
ity of a Feshbach resonance (see Fig. |]) . Two parameters 
characterize the Feshbach resonance: (i) the energy mis- 
match e between the bound state and the zero energy 
edge of the continuum states for the colliding atom pair, 
and (ii) the inverse lifetime k of the bound state. As 
an alternative, photoassociation may be used to directly 
generate the resonant coupling [jll 12 1. In photoassocia- 
tion, a two-photon Raman transition is used to couple the 
atomic continuum states and a specific bound molecular 



level. In that case e represents the detuning energy of the 
Raman lasers from the atom-molecule transition, and k 
denotes the two-photon Rabi frequency proportional to 
the laser intensities and to the usual overlap integrals. 



Closed Channel 




FIG. 1. Feshbach resonance. Atoms collide with rela- 
tive kinetic energy near zero as indicated by the dotted line. 
This energy is quasi-degenerate with a bound state in a closed 
channel potential. Typically, e can be tuned experimentally 
by changing the magnitude of an external magnetic field. 

An effective Hamiltonian for this coupled atom- 
molecule system may be written as 

H = J d 3 xtPl{x)H a (x)^ a (x) + ^ m {x)H m {x)vb m {x) 



+ 



d 6 X 4>l {x)lpl (x) 1p a (x) 1p a (x) 



+ U am I d 3 x^l(x)^l l (x)rl> m (x)ij> a (x) 



d 3 X #h{ x wL( x Wm(&Wrn(&) 



d 3 xtpj n (x)ip a (x)^ a (x) + h.c. , 



(1) 



where ij) a [x) and ijj m (x) are bosonic field operators which 
annihilate an atom or molecule respectively at coordi- 
nate x. Atomic binary collisions are characterized by 
the mean- field energy per unit density, ?7 aa = Airh 2 a/m, 
where m is the atomic mass, and a is the scattering length 
for atom-atom collisions. Equivalent definitions apply for 
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C/ am and U mm corresponding to atom-molecule collisions 
and molecule-molecule collisions. The process of inter- 
conversion of atom pairs into molecules is characterized 
by g — V K U aa M. The free Hamiltonians are 



and coupling elements 



H a (x) = - — Vl + V a (x) - fl a 
2m 



H m {x) 



4m 



Vl + V m (x)- fi m , 



(2) 



where we include V a (x) and V m (x) for generality to allow 
for the possibility of external potentials. The chemical 
potentials of the atomic and molecular fields are denoted 
by \i a and /i TO respectively, such that e = fi m — 2fi a . 

We solve for the evolution of the fields, assuming the 
principle of attenuation of many-particle correlations in 
the dilute gas. That is, we find the expectation value 
of the Heisenberg equations for ip a (x) and ip m (x) and 
separate explicitly the mean values of the field operators 

(t>a(x) = (lpa{x)) 

4>rn{x) = (tp,n{x)) (3) 

from their zero-mean fluctuating components 



Xa(x) = ljja(x) - tf> a (x) 
Xrn(x) = lp m (x) - (f> m (x). 



(4) 



In order to close the dynamic equations, we expand prod- 
ucts of three of more fluctuating operators using Wick's 
theorem (i.e. we assume a Gaussian reference distribu- 
tion for the many-body quantum state). Since we intend 
to determine the pairing statistics of atoms and do not 
anticipate intrinsic pairing of molecules to play a signif- 
icant role, we assume a classical field for the molecular 
quantum state (i.e. we specify the molecular field as an 
exact eigenstate of ip m (x)). Following this procedure, we 
derive a description of the coupled atom-molecule sys- 
tem, which involves the following dynamical quantities: 
the condensates (j) a {x) and 4> m (x), and the atomic fluc- 
tuations as described by normal densities GN(x,y), and 
anomalous densities G A (x,y), defined by 



G N (x,y) = (xi(y)Xa(x)) 
G A (x,y) = (Xa(y)Xa(x)) . 



(5) 



Following this prescription, we arrive at the following dy- 
namical equations for the condensates: 

ih ^a(x)_ = ( Ha{x) + , J , a (x))^ a (x) + e( X )p a (x) 



ih 



dt 
d(f) m (x) 
dt 



(H m (x) + r m (x))(t> m (x)+$(x), (6) 



with corresponding mean-field potentials 

t a (x) = U aa .(\(/) a (x)\ 2 + 2G N (x,x)) + U m \4> 

t m (x) = U &m (\4> a {x)\ 2 + G N {x,x)) + U mm \4> m {x)\ 2 , (7) 



0(as) = U aa G A (x,x) + g4> m {x) 
$(x) = |(0 a (x) 2 + G A (x,x)). 



(8) 



A concise representation for the evolution of the atomic 
fluctuations is given by 



= ^g- est , 

dt 

which is written in terms of the matrix 



(9) 



y{ ' yj " 1 (xi(y)xUx)) (Uy)xl(x)) J ' (1Uj 



and Bogoliubov self-energy of general structure 



E(x,») = 



H{x,y) 0{x,y) 
-Q(x,y)* -H(x, V y 



(11) 



The 2x2 elements of E(a;, y) are each diagonal in x and 
y due to the assumption of a contact potential for all 
scattering diagrams in Eq. (|]J). Accordingly, by defining 
H(x,x) = H a (x) + T(x,x), the full solution is 

f(x, x) = 2U^(\4> a (x)\ 2 + G N {x, x)) + U &m \^ m (x)\ 2 
Q(x,x) = U^aix) 2 +G A (x,x)) +g<(> m (x). (12) 

This forms a universal description of the Hartree-Fock- 
Bogoliubov theory of the coupled atom-molecule conden- 
sate system. As an example of the implementation of 
these equations, we now apply the theory to the case 
of a uniform gas with the assumption that we may ig- 
nore the mean-field potentials (i.e. we set to zero C/ aa , 
C/ am , and U min ), along with the external potentials V a {x) 
and V m (x). Due to translational symmetry, the result- 
ing atomic and molecular condensates are spatially ho- 
mogeneous, and the normal and anomalous densities, 
Giq{x,y) and GA{x,y), depend only on the magnitude 
of the relative coordinate, i.e. \x — y\. 

We convert the equations to dimensionless form in 
the following manner. We define the number density 
n as the number of atoms plus twice the number of 
molecules per unit volume. Then g^/n denotes a char- 
acteristic coupling energy of the atomic and molecular 
mean-fields. A dimensionless time r can therefore be 
defined by r = g^/nt/h. A dimensionless coordinate 
r = \x — y|/2£ is associated with the length scale ( corre- 
sponding to the formation of molecules from atom pairs. 
This is effectively a "healing length" for molecule forma- 
tion found by balancing the relative kinetic energy with 
the resonance coupling energy, i.e. h 2 /2/i£ 2 = g^/n where 
[i = m/2 is the reduced mass. Making systematic sub- 
stitutions 4> a = 4> a (x) / y/n, tfirn = <j> m {x)/ y/n, G N (r) = 

G N (x,y)/n, G A {r) = G A {x,y)/n, and A = e/gy/n gives 
the complete system of equations in dimensionless form 
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dr 



dT 

dG N (r) 

1 dT 

dG A {r) 



A + \{<t>l + G A {r)\ r = ) 



b m G A (r)-r m G A (r) 

-V 2 r G A {r)+^ n (2G N {r)+n5^{r)) , (13) 

where 5^ (r) is the isotropic three-dimensional Dirac 
delta function and we have defined 77 = as the 

inverse diluteness parameter of the coupled atomic and 
molecular gas. There are explicitly conserved quantities 
in these equations corresponding to normalization (den- 
sity of atoms plus twice the density of molecules) i.e. 



\<j>a\ 2 + G N (r)\ r=Q 



= 1. 



(14) 



and energy density u (energy due to detuning and cou- 
pling plus the kinetic energy of the normal gas) 



u = -A\<j>„ 



Re[C(^ + G4r)| r=0 )] 



"2 w)U 



(15) 



Interestingly, Eqns. (|T^) cannot be directly integrated 
as written. The delta function describes the spontaneous 
breakup of molecules into atom pairs. This process is 
reminiscent of the decay of an excited atom leading to 
the spontaneous emission of a photon. In that case, the 
interaction with the unoccupied vacuum modes of the ra- 
diation field gives both an energy width and an energy 
shift (Lamb shift) to the decaying state. Analogously 
in the atom-molecule system, there is an energy shift of 
the molecular state containing contributions from all dia- 
grams representing the breakup of molecules into virtual 
atom pairs. The virtual atom pairs may form an inter- 
mediate state off the energy shell. Summing over the full 
spectrum of wave numbers for the pair gives an infinite 
shift to the molecular level. 

This divergence is reconciled by noting that we have 
been inconsistent in not including the self-energy of the 
molecular quantum state in the original definition of the 
chemical potential fi m . We carry out the renormalization 
in the following manner. We place an artificial bound on 
the momenta of the atom pair by replacing the delta 
function by a three dimensional Gaussian of standard 
width ay normalized to have unit volume 



(16) 



Physically, this accounts for the fact that the real poten- 
tial is not exactly a contact potential. The resulting en- 
ergy shift is then finite and we derive its value according 
to the following prescription. We define a three dimen- 
sional Fourier transform 



G A {k) = / d 3 rG A (r)e 



-ir-k 



(17) 



which for the isotropic case is simplified to 

G A {k) = ^- rG A (r)sm(kr)dr. (18) 
* Jo 

For sufficiently high wave numbers, an approximate equa- 
tion may then be written for G A {k) using Eq. (pj 



i^^- = k 2 G A (k) + r ] (k)^ m , 
(It 



(19) 



where rj(k) — nexp(—a 2 k 2 /2). Taking G A (k)\ T =o = 0, 
the solution of this equation is 



G A {k) = -i V (k) / <MT')e lfe2(T '~ T) dT' 
Jo 

~ "■p 7 ?( fc )^»i( T ) • 



(20) 



Calculating the energy shift of the molecular level re- 
quires substituting this result into the evolution equa- 
tion for d<p m /dT in Eq. (|l3|). Using the Fourier integral 
G A (r)| r = = (27T 2 )- 1 f™k 2 G A (k)dk, the resulting shift 
may be incorporated into a renormalized detuning A 



A- 



77 V 2ir 



(21) 



The end result of this analysis is that Eqns. (|l^) are satis- 
factorily renormalized by making the two substitutions, 
Eq. (|lG|) and Eq. (|2l|). The evolution of any relevant 
observable will then be independent of the choice of ay, 
providing ay is chosen to be sufficiently small. 

As a numerical example, we examine the 23 Na Fes- 
hbach resonance at 907 G ||. For this system k = 
2irh x 44 MHz, and we consider the bound state to be 
on-resonance with the colliding atom pairs. Other pa- 
rameters are the scattering length a — 60 do , where ao is 
the Bohr radius, and we take a density of 10 15 cm -3 . This 
gives an inverse diluteness parameter of r\ = 16.5. In the 
simulation in Fig. |^ we show the evolution of the atomic 
and molecular condensate fields, comparing the Hartree- 
Fock-Bogoliubov theory we have derived with the pre- 
dictions of mean-field theory. The initial condition is 
taken as 95% of the population in the atomic conden- 
sate and 5% in the molecular condensate. The persistent 
large scale oscillations of the population of the atomic 
and molecular condensates as seen in a solely mean-field 
theory dampen out when the coupling to the normal gas 
is included. This is partly due to the fact that when 
molecules are formed from atom pairs, the pair correla- 
tion function at that point is depleted, and partly due to 
the possibility for spontaneous breakup of the generated 
molecules into the normal component of the gas. 

In Fig. |^, we illustrate the behavior of the normal com- 
ponent during this simulation. The density of the normal 
gas is Gjv(r)| r= o. The temperature of the normal gas is 
found from the spatial scale over which this correlation 
function decays. For reference, for an equilibrium classi- 
cal gas, Gat is Gaussian with a standard deviation given 
by the thermal de Broglie wavelength. 
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FIG. 2. Time evolution of the atomic condensate \(f> a \ 2 
and molecular condensate 2\4> m \ 2 for solely mean- field theory 
(dashed) and for the Hartree-Fock-Bogoliubov theory (solid). 




FIG. 3. Evolution of the normal density. 




FIG. 4. Evolution of the anomalous density. 



In Fig. |4| we show the evolution of the anomalous 
fluctuations for this example. For a classical gas, the 
anomalous fluctuations are explicitly zero, so clearly the 
quantum statistics of the normal component generated 
here are not the usual classical thermodynamic equi- 
librium. One of the reasons for this is that there can 
never be an odd number of atoms in the thermal cloud 
since atoms are spontaneously generated from molecules 
in pairs. This situation is similar to the formation of a 
squeezed vacuum in optics using a laser pump to drive 
a parametric amplifier in order to produce photon pairs. 
Note the value of |Gyi(r)| near r — depends on a" 1 and 
is not an observable. 

In conclusion, the damping of the coherent atom- 
molecule oscillations is accompanied by an increase in 
density of the normal component. The behavior shown 
in Fig. H resolves the conceptual difficulty associated with 
the rapid atom-molecule oscillations predicted by mean- 
field theory — how can pairs of atoms find each other 
rapidly enough to form molecules at the predicted rate? 
The proper treatment of correlations shows that the os- 
cillations are slower in frequency and moreover rapidly 
damped. 
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